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Pinning-induced pn junction formation in low-bandgap two-dimensional 
semiconducting systems 
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A model is presented for pn junction formation near metal-semiconductor contacts in 
two-dimensional semiconducting systems such as graphene. Carrier type switching 
occurs in a region near the metal-semiconductor junction when energy band bending 
leads to a crossing between the junction Fermi level and the Dirac energy. A bias- 
dependent depletion region occurs due to the minimization of carrier density, which 
is shown to act as an additional parasitic resistance in devices. The pn junction 
resistance is demonstrated by its implementation in a transfer length structure. 
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The advent of two-dimensional (2D) semiconducting systems, such as graphene and the 
transition-metal dichalcogenides, has driven significant efforts to utilize these systems for 
electronic devices.^ In the vast majority of embodiments, electronic devices like the field 
effect transistor and various sensor designs utilize a modulated channelP In those cases, 
metallic electrodes are used to make ohmic contact to the semiconducting channel. Most 
studies address the minimization of parasitic access resistance via the contact resistance {Rc) 
of the contacts to graphene.®® However, little attention has been given to the manifestation 
of spuriously high-resistance pn junctions formed in close proximity to the contact. The 
formation of the pn junction is due to the contrast between the Fermi-energy pinning at 
the metal-semiconductor interface and opposite carrier typing outside the contact region 
that results from ambient electrostatic bias. Ultimately, the pn junction resistance results 
as an additive parasitic resistance to the other parasitic resistances (access. Race-, and Rc) 
present, which may subsequently serve to impair device performance. This work presents the 
development of an electrostatic model that describes the formation of such a pn junction 
with low carrier density in proximity to a metal-2D semiconductor contact. The energy 
intersection between the pinned Fermi level under the metal-semiconductor contact and the 
2D semiconductor Dirac energy {Ed) is central to the development and results in the low 
carrier density bias-dependent pn junction. The model includes calculations for the width 
of minimum conductivity, its dependence on voltage bias, and how it is manifested in a 
conventional test structure. More complex effects such as Klein tunneling and the role of 
band chirality on carrier generation-recombination are not incorporated in the electrostatic 
model but likely play a role in the transport picture.® 

Graphene and some transition-metal dichalcogenide semiconductors exhibit a unique lin¬ 
ear energy dispersion relation due to hybridized sp^ covalent bonding within their bravais 
lattice. In this work zero-bandgap graphene is used as an example from the broad 2D 
semiconducting families. Figure shows an illustration of a graphene device cross-section 
near one ohmic contact (a) with the linearized energy band diagram (b) and resultant total 
charge distribution (c) for several voltage biases. The ohmic contact was assumed to pin 
the graphene Fermi level at an energy of Pp only under the contact.®^ The linearization 
of the energy band diagram is an approximation made only in a small energy range about 
the Fermi level-Dirac energy (F-D) intersection where total charge density is minimized and 
resistance has a maximum. Realistically, the Dirac energy (or conduction and valence band 
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intersection) bends throughout the length of the device in the x-direction in coherence with 
a change in charge density (Poisson). However, for the small energy range about the F-D 
crossing the linearized energy is sufficient to describe the carrier typing and therefore, the 
electrostatics of the system in the vicinity of the ohmic contact. 

An approximation has been made for the resistance associated with the region beyond 
the linear energy region where |(i?D — Ef)\ 3> \rip\. In this region, the total charge density 
exponentially increases according to Fermi-Dirac statistics and yields a resistance much less 
than in the region where the F-D crossing occurs. Therefore, this low-resistance extension 
has been treated as the sheet resistance outside the proximity of the F-D intersection (regions 
beyond in Fig. 1 (c)). The pinning energy polarity at the metal-graphene junction and 
the ambient electrostatic bias away from the contact are responsible for setting up the 
conditions leading to the F-D intersection and therefore, the attributes of the pn junction. 
Those conditions determined the polarity of the Ed slope, number of F-D intersection within 
a given region, and length of total charge minimum referred to as the ‘depletion width’ {Wd) 
in this work. If two F-D intersections occur, the total resistance becomes approximately 
twice the resistance of a single pn junction. 

The carrier distribution in graphene results from the linear dispersion relation and there¬ 
fore follows from the full Fermi-Dirac integral, n,p{p) = qNGEi{±p) where Ei{i:p) = 
(F(z -I- 1))“^ M*(l -f The longitudinally-dependent total charge is given 
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where Nq = {gsgvl‘^T^){kBT/hupY is the 2D density of states with gg and g^ are the spin 
and valley degeneracies, respectively, vp is the Fermi velocity, q the elemental charge, F(...) 
is the gamma function, h is the reduced Planck’s constant, ((>p{x) = ksT^p^x) is the energy 
along the channel normalized by the thermal energy kpT, kp is Boltzmann’s constant, and T 
the temperature. The linearized energy band was approximated as, E — Ef = Va{x — Xo)/L, 
where Xo is at the ohmic contact edge as illustrated in Fig. The normalized longitudinal 
energy term used in Eq. [^follows from the boundary conditions of the linear energy band. 
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where L is the length of the graphene sheet between the two biased contacts as dehned in 
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FIG. 1. (a) Illustration of the 2D structure, (b) lateral band diagram in graphene along the source- 
drain direction at several different voltage biases, and (c) the resultant carrier density profiles along 
the graphene channel. 
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Fig. i4 is the applied voltage, and rjp is the contact pinning energy. The depletion width 
is dehned as twice the distance from contact pinning to the Fermi level crossing and follows 

qVa 

In the conventional case of a three-dimensional semicondnctor pn jnnction, the depletion 
width is dehned by the balance of the internal held and the dihusion of free carriers. The 
linear energy approximation acconnts for the complete electrostatic pictnre for the pinning- 
indnced pn junction in graphene. The total resistance of the depletion width accounting for 
the single Fermi level crossing, is given by the integral relationship. 


Wa Wa 

^ 1 r dx ^ f 

Wg J a{x) qWg J pnn{x) + 
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where a{x) = q{pnn{x) -|- ppp{x)) is the graphene conductivity and is dependent on the areal 
charge densities, n and p, and the respective carrier mobilities, pn and Pp, and Wg is the 
cross-sectional width of the graphene sheet. 

Parametric calculations are shown in Fig. [^wherein Figs, [^(a) and (c) show total charge 
density and correspondent resistivity {p{x) = 1/ct(x)) in the graphene sheet, respectively, 
as a function of distance away from the ohmic contact for several longitudinal voltages. A 
pinning energy of 4 eV has been used, which is commensurate with values reported.®® The 
charge density and resistivity minima progressively decreased with an increase in applied 
voltage as the F-D intersection migrated toward the ohmic contact. The rapidly diminished 
value of p away from the F-D crossing (total charge minima) validates the previously stated 
approximation. Figure (b) and (d) show depletion width and depletion region resistance 
in the graphene sheet, respectively, as a function of applied voltage for several values of 
pinning energy. As a greater bias voltage is applied, energy-band bending is enhanced and 
the slope of the linear region becomes greater leading to the trends shown in Fig. The 
bias-dependence of Rd is of note within the context that parasitic resistances are typically 
hxed (e.g. or Race)- A dynamic parasitic resistance like Rd may serve to complicate 
the extraction of other intrinsic device metrics (i.e. transistor transconductance, frequency 
performance, etc.). One method that may prove to disable the formation of the pn junction 
is doping graphene. By doping, the Fermi level in graphene can be signihcantly shifted such 
that no F-D intersection may occur. Both Wd and Rd linearly decreased with increasing 
voltage following Eq. and respectively. 
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FIG. 2. Parametric characteristics for the pn junctions formed in graphene in proximity to an 
ohmic contact; (a) Charge density versus longitudinal distance, (b) resistivity versus distance for 
various applied voltages, (c) depletion width versus applied longitudinal voltage, and (d) depletion 
region resistance versus applied voltage for various pinning values. 

The effect of the pn junction on a device is illustrated through implementation into 
a transfer length structure (TLM). The TLM structure is a conventional test structure to 
extract information on R^. and Rgh- The total resistance of a TLM structure is conventionally 
given by R{L) = 2Rf. + ^Rsh, where Rgh is the sheet resistance. For the scenario considered 
here, we assume there is a single pn junction between TLM contacts and the depletion width 
is small compared to the separation between contacts, Wd -C L. The latter assumption is 
verified by the calculation of Wd as shown in Fig. [^b). These assumptions result in the 
TLM resistance relationship for graphene, 

RtLM = 2Rc + Rd + Rsh ^ (^) 

with Rd being the depletion region resistance as given by Eq. |^and a is the distance between 
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FIG. 3. Resistance of TLM structure comprised by graphene with and without the inclusion of 
a single pn junction showing a corresponding resistance increase. The inset shows the graphene 
TLM configuration and ID linearized energy band diagram between two contacts. 

the contacts and the depletion width as illustrated in the inset in Fig. Written in terms 
of the longitudinal graphene conductivity and invoking the constraint, L S> a, the total 
resistance of a graphene TLM structure with a single pn junction follows, 

Wa 

T ^ r 

R{L) = 2R^ + —R,, + W ^6) 

0 

Rd 

Equation shows that the presence of a single pn junction in the graphene sheet adds a 
resistance proportional to the total bias-dependent resistivity of the junction to the TLM 
resistance. The pn junction resistance is independent of separation, L, and therefore is 
additive. The manifestation of additional pn junction resistance is demonstrated in Fig. 
where a single pn junction is included in a TLM resistance curve with a Fermi pinning of 4 eV 
compared to a curve that does not include a pinning-induced pn junction. Other pertinent 
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graphene parameters used in the calculation are Wg = 20 pm, Rsh = 300 hl/D, i^n = = 

1000 cm?jVs, T = 300K, and vp = 1.1 x 10® cmjs. The linear energy approximation should 
hold for other 2D low bandgap semiconductors. 

An electrostatic model has been presented for pn junction formation in graphene in close 
proximity to ohmic contacts where Fermi level pinning occurs. It has been shown that 
electrostatic depletion of charge within the pn junction introduces a bias-dependent high- 
resistance region. The effect of pinning-induced pn junction formation in graphene has been 
illustrated for the case of a TLM structure, wherein the junction resistance was shown to be 
an additive term to the total TLM resistance. The broader impact of pn junction formation 
in low-bandgap two-dimensional systems may be spurious resistance that impairs device 
performance. 

We would like to thank J. G. Champlain at the Naval Research Laboratory for helpful 
technical discussion. 
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